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I.  Planning  and  Survey  Stage 


The  work  began  with  a  computer  search  for  relevant  material  published 
during  the  last  two  years  in  the  Chemical  Abstracts  and  Physics  Abstracts.  From 
the  Chemical  Abstracts  a  total  of  113  references  responded  to  a  call  from  the  key 
words  "electromigration"  and  "electrotransport".  These  were  surveyed  for 
applicability  and  availability  and  26  rated  high  on  both  counts.  The  corresponding 
reference  list  was  much  longer  for  the  Physics  Abstracts.  It  was  decided  to 
require  also  mention  of  "thin  film"  as  a  key  phrase  and  this  reduced  the  list 
substantially  -  to  25.  All  seemed  pertinent  and  were  provided  with  substantial 
abstracts.  There  was  naturally  a  considerable  degree  of  overlap  between  the  two 
lists.  This  background  has  proved  most  useful  in  selecting  a  topic  for  in  depth 
study  and  also  in  giving  a  general  picture  of  the  most  recent  aspect  of 
electromigration  problems. 

II.  Generation  of  Grain  Boundary  Network. 

In  the  course  of  this  literature  survey  there  surfaced  one  particularly 
pertinent  paper  by  K.  Nikawa*,  which  described  a  detailed  computer  simulation, 
using  Monte  Carlo  methods,  of  a  rather  simplified  model  of  stripe  deterioration  by 
grain  boundary  transport.  The  accompanying  reference  list  ^_5  includes  several 
earlier  papers  in  the  general  field.  It  was  in  this  area  we  decided  to  concentrate 
our  computing  efforts. 

The  first  step  in  this  procedure  was  to  devise  a  system  for  generating  a 
realistic  grain  boundary  network.  This  was  planned  as  a  two  step  process. 

A.  First  a  random  network  of  grains  was  generated  according  to  the 

Voronoi  prescription,  which  has  been  discussed  in  a  review  article  by  Weaire  and 
6 

Rivier  .  The  method  is  to  establish  a  random  placing  of  points  in  two  dimensions. 
Next  one  constructs  the  perpendicular  bisectors  between  all  pairs  of  adjacent 
points  and  portions  of  these  lines  now  constitute  the  grain  boundaries  of  the 
network.  The  physical  interpretation  could  be  that  the  random  points  correspond  to 
nucleation  centers  for  crystallization  which  grow  uniformly  until  they  meet  at  the 


perpendicular  bisectors,  i.e.  grain  boundaries.  In  programming  this  concept  for  a 

7 

stripe  of  a  particular  width  we  made  use  of  a  prior  program  by  I.  K.  Crain  for  a 
limited  network  and  were  successful  in  generating  reasonable  structures  for 
different  size  grains.  As  a  last  step  grain  orientations  were  established  by  Monte 
Carlo  selection. 

B.  Next  an  effort  was  made  to  equilibrate  the  intergranular  tensions  at 

the  vertices  of  the  network.  We  assumed  that  these  tensions  would  practically 

equal  to  the  intergranular  energies  (on  the  assumption  of  negligible  temperature 

dependence).  For  these  energies  a  dislocation  model  was  again  invoked  which  set  the 
grain  boundary  energy  proportional  to  the  dislocation  density  or  sin  §  for  6  <  37° 

and  kept  the  energy  constant  for  9  >  37°,  where  9  is  the  angle  of  misorientation 

between  the  adjacent  grains.  Since  the  grain  boundary  diffusivity  also  varies  in 

somewhat  the  same  way  as  the  tension,  it  seemed  desirable,  perhaps  even  important, 

to  try  to  establish  tension  equilibration.  This  proved,  however,  quite  a  difficult 

business,  as  might  have  been  anticipated  for  a  multidimensional  search  for  a 

equilibrium  configuration.  For  a  time  this  effort  was  shelved  in  favor  of  an 

exploration  of  the  actual  mass  motions  with  the  grains. 

III.  Modelling  the  Electromigration  Process 

A.  Characteristics  of  the  grain  boundaries 

The  key  quantity  in  the  simulation  process  is  m0(x,t),  the  amount  of 
matter  in  grain  boundary  segment  (denoted  by  subscript  a)  per  unit  length. 

Initially  all  grain  boundaries  are  assumed  ungrooved  and  the  m0  have  a  constant 
value  biq  *  Nd06  where  N  is  the  number  of  atoms  per  unit  volume,  dD  is  the  thickness 
of  the  stripe  film  and  6  is  the  (nominal)  width  of  the  grain  boundary  segment.  In 
some  regions  m0  >  idq  and  there  is  matter  build-up  which  continues  until  iHq  ~  m^  at 
which  point  extrusion  starts.  In  contrast  to  those  grain  boundaries  which  show 
matter  excess  there  are  others  which  show  matter  deficit  m0  <  b^.  Because  these 
regions  are  critical  in  void  formation  their  treatment  is  more  detailed  and 
complex.  Matter  deficit  in  the  grain  boundaries  (vacancies)  tends  to  condense 


as  grain  boundary  grooves  at  the  surface.  As  we  show  later,  the  number  of 
vacancies  still  in  the  grain  boundary  is  relatively  small  after  a  short  initial 

period.  The  matter  missing  in  the  grain  boundary  grooves,  per  unit  length  of 

2 

grain  boundary,  is  taken  to  be  Ndg  tan  3  (tan  6  is  small  and  has  been  arbitrarily 
set  at  0.1)  where  dg  is  the  depth  of  the  groove  and  3  is  the  angle  at  the  groove 

base. 


With  the  appearance  of  the  grooves  another  mechanism  for  mass  transport 

becomes  operative  -  surface  motion  down  the  groove  walls.  (Presumably  surface 

diffusion  over  the  upper  surface  of  the  stripe  is  inhibited  by  a  passivating  layer.) 

This  process  should  proceed  rapidly  and  be  independent  of  the  original  grain 

misorientation  at  the  boundary.  Its  effect  is  to  widen  the  groove  and  therefore  a 

new  quantity,  w(x,T),  is  introduced  which  denotes  the  grain  boundary  width  at  the 

surface  for  each  point  along  the  segment.  This  action  greatly  increases  the 

effective  3  so  that  3  >  3'  where  tan  3*  “  — ^ — 

g 

In  addition  to  the  m(x,t)  quantities  to  describe  the  state  of  the 
segment,  one  might  also  consider  g(x,t),  the  groove  volume  per  unit  length.  It 
appears,  however,  that,  in  general,  g(x,t)  *  (m 0  -  m(x,t)  will  be  close  to  unity, 
since  the  vacancy  volume  in  the  grain  boundary  is  relatively  small. 

B.  Equations  for  mass  motion 


Within  the  grain  boundary  the  appropriate  mass  transport  equation  is 
3m(x,t)  3  ,  r  3m(x,t )  .  ,  .  , i , 

at  -  *T  {D0i  t  3x  -  A  B(x»t)J)  (1) 


Distance  along  the  grain  boundary  segment  is  denoted  by  x.  The  grain  boundary 
diffisivity  is  D0^  where  a  again  numbers  the  segment.  We  show  in  next  section 
D^0(9)^  where  9  is  the  misorientation  angle  between  the  grains  meeting  to  form  the 
boundary.  The  coefficient  of  the  driving  force  term,  A,  is  Aq  cos  4*  where  $  is 

the  angle  between  the  grain  boundary  and  direction  of  E  and  Aq  is  J ^ —  about 

O.lym  for  E  *  2.2V/cm,  Z*  ■  20  and  T  "  150°  C.  All  changes  in  m  are  treated  on  an 


equal  basis,  whether  they  arise  from  changes  of  vacancy  concentration  in  the 


grain  boundary  or  a  gradient  in  groove  depth. 

Preparatory  to  setting  the  equation  up  for  computer  application  each 
segment  had  its  length  divided  into  a  discrete  number  of  evenly  spaced  check 
points.  The  spacing  is  denoted  by  a^  and  is  roughly  about  0.2ym  for  all  segments. 
Eq  (1)  is  now  written 

3m(x,t)  =  °lo  |-mn+l  -  ^mn  +  mn-l  -  mn+l  -  mn-l  .  i  /,  \ 

at  —  l - ; -  - 2 -  V  <u) 


Where  the  digitized  values  fo* 


have  been  used  on  the  right.  At  the  ends  of 


the  segments  (n  =  0  or  nt]  ,  equation  is  slightly  altered. 


3m(0,t)  =  1  IQ  rml 

3t  a  a  *■ 


A  m,  ] 
o  1 1 


In  a  similar  way  the  surface  motion  on  the  groove  walls  can  be  expressed 
in  terms  of  g(x,  t)  *  N  x  volume  of  the  groove /length. 

8a^Lii'-r-  Ds  gu.t>]  (2) 

t  x  x 

Where  Ds  is  the  diffusion  on  the  free  surface  and  Ag  is  the  surface  driving  force. 


As  before  the  digitized  version  is 

1&.  J9a  rSntl £gn— t-J 

3t  “  a  1  a 


i—l  —  Ac 


and  for  the  segment  ends 

ifia  _  ^  [SI  ~  & a_  Aggl]  (2b) 

3 1  a  a 

(In  a  later  section  we  give  an  argument  why  the  diffusion  term  should  be  omitted 
from  eq.  2.) 

C.  Assumed  stripe  structure 

It  is  reported  that  most  metallized  stripes,  particularly  for 
aluminum,  exhibit  a  pronounced  (111)  structure.  Such  a  structure  implies  that 


pure  tilt  interfaces  constitute  most  of  the  grain  boundaries.  If  we  invoke  a 
dislocation  model  for  the  small  angle  grain  boundaries,  these  would  then  consist  of 
edge  dislocations  directed  perpendicular  to  the  substrate,  hence  the  notation  D^0 
for  diffusion  in  the  stripe,  perpendicular  to  the  dislocations.  Accordingly  the 
assumption  is  made  that  for  large  angle  grain  boundaries  (0>37°  as  a  nominal 
dividing  line),  diffusion  is  isotropic  at  a  value  Dj.  Holding  to  the  dislocation 
model  for  grain  boundaries  one  obtains  for  the  diffusion  parallel  to  the 
dislocation  and  directed  perpendicular  to  the  substrate 

Dj  *  D^  sin  0/2/sin  37°/2  for  0*  37°  (3) 


For  diffusion  along  the  stripe  one  modifies  a  formula  of  Li  to  get 


(4) 


where  D0  is  the  diffusion  through  the  bulk  and  a0  «  rD/h.  Here  rQ  is  the 

dislocation  core  radius,  b/h  =  2  sin  ^  and  b  is  the  Burgers  vector.  (The  original 

2 


1  = 

Li  formula  was  —  — - 

U  l)q 


(l-a0)Do  +  OqDi 


).  Nominally  a 


sin  0/2 
sin  37°/2  . 


D.  Difficulty  with  Diffusion 

The  direct  integration  of  equs  (1)  and  (2)  poses  a  problem  in  that 
useful  convergence  of  the  diffusion  part  requires  some  sublety  in  method  and  a 
very  considerable  number  of  time  steps.  In  terms  of  an  average  segment  length  L^v 

the  natural  time  constants  of  the  system  are  [d^  L~2  ]  (about  .6  x  l(Ps)  and 

[Dj  A  Lav“1 ]  (about  5  x  lO^s).  As  both  of  these  constants  are  very  short,  of  the 
order  of  an  hour,  it  is  clear  that  the  approach  to  void  formation  must  not  be 
limited  to  time  steps  of  such  small  size.  Accordingly  a  procedure  was  developed  to 


integrate  over  a  longer  period. 
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E.  Steps  in  the  program  simulation 


The  individual  stages  that  are  involved  in  each  time  step  are  here 


described  in  order,  somewhat  after  procedure  used  in  the  twelfth  status  report. 


a.)  First  stage  -  matter  accumulation  (or  depletion)  at  each  vertex, 


Fv,  in  accord  with  eq  (lb). 


,m  -m  . 

D  [  -Am,], 

o  ol  L  a  a  1 J 


where  the  o-summation  is  over  the  segments  which  meet  at  the  vertex  in  question. 


For  the  first  time-step  all  the  m's  are  equal  and  only  the  A0  term  contributes. 


b.)  Redistributing  Fv  among  the  segments  meeting  at  V. 


For  each  segment  a  portion  of  F,  namely  f0,  is  redistributed  as  a 


flux  coming  into  the  segment. 


f  =  F  D 


v  oj. 


/  Z  D, 


c.)  Method  of  redistribution 


The  thin  film  solution  for  the  diffusion  equation  is  used,  working 


from  both  ends  of  the  segment, 
3m(xjt )  _  , 


-x  /4D  t 
e  la 


(it  D  t) 
la 


One  would  like  to  be  able  to  state  what  m(x,t)  is  at  time  t. 


Integration  gives  m(x,t^)  =  (ir)  f^  x  5/ 


j  Where  £(y)  =  y  -  /—  erfc  y* 


This  integration  technique  allows  one  to  greatly  increase  the  time,  from  lO^  to  10^ 


seconds.  For  such  long  times,  however,  the  diffusion  thin-film  solution  will 


usually  have  a  large  value  at  the  opposite  end  of  the  segment.  These  "tails"  have 


to  be  evaluated  by  integrating  m(x)  from  xj  to  «, 


/  m(x)  dx 


t  n(y  )  *  T(x  ,t  ) 
111  11 


■\V V 


:-v 


k 


v  .N  .  .«  •• 


One  then  sums 


where  n(y,  )  =  (1  +  2y  )  /  00  e  y  dy  -  y,  e^l  and  y,  =  2  /  D  .t 
1  1  y !  1  1  O1  1 

the  KL^.t)  at  each  vertex  for  the  "tail”  contribution  from  all  the  segments  that 

meet  there.  This  gives  a  new  which  is  then  iterated  to  give  smaller 

contributions  to  the  m0(x,ti)  the  network. 

d.)  After  these  "tailing"  corrections  have  converged,  the  next  stage  is 

to  consider  the  grooving  that  occurs  near  those  vertices  where  a  material  deficit 

has  developed.  If  m(x,t, )  -  m  <  0,  it  is  reasonable  to  expect  large  fraction  of 

1  o 

this  deficit  will  be  converted  into  the  grain  boundary  groove  volume,  g(x,  t^),  or 


g (x , t ^ )  =  K(x,t1)  {mQ  -  m(x, 1 1 ) }  =  K(x,ti)  mi(x,ti)  (8) 

K(x,tj[)  can  be  estimated  from  a  model  wherein  the  diffusion  of  vacancies 
in  the  z-direction  (upward  from  the  substrate)  is  explored.  Suppose  m(x,z,t^)  * 

m(x,tj)  Z(z)  where 


Z(z) 


cos 

2db  2db 


(9) 


Here  z  =  0  at  the  substrate  and  z  *  d^  at  the  bottom  of  the  groove,  d^  =  dQ-  dg 
where  dg  is  the  groove  depth.  Introduce  now  T(x,ti),  t-^ie  number  of  vacancies  still 
in  the  grain  boundary,  T(x,ti)  =  mi(x,t})  ~  g(x,ti) 


3 


T(x,t) 

3t 


3m  (x,t) 

-A^T(x,t)  + 


(10) 


TT 


i 


Here  A  =  D a ( 0 ) 

substrate  and  ■ 


( 2d  )  from  the  equation 
b 


3m^ 

IT 


could  be  taken  from 


for  diffusion  perpendicular  to  the 


(eq  6)  -  if  tailing  correction  is  neglected. 


To  get  a  rough  solution  disregard  the  Gaussian  factor  since  long  times  are 


3m  1 


-1 

2 


1  J_ 

2  2 , 


involved.  Roughly  then  ^  “At  and  m  ~  2A(t^  -  t^  )  when  tm  is  a  short 


time  cut-off  which  replaces  the  exponential. 


3T  IT  +  4  2 

3T=  'AT  +  Ae 


(10a) 


Discarding  the  homogeneous  solution  ce  we  have  for  T 


-At  t  1  t  1 

T(t)  =  Ae  l!  1  et>X  (ti-t')~2  dt  ’  =  A  /  1  e~Xr  r~2  dr 
n 


T(t)  -  2A  x  [erf (X  Z  ti>  -  erf(X^  tm) ]  “  2A  X 
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and  finally  K(x,t)  =  (1  -  (Xt^)  }  (11) 

Often  K(t)  is  close  enough  to  unity  to  neglect  the  difference.  Actually  the  form 


3m, 


of  (9)  is  quite  independent  of  the  formula  for  .  To  determine  dg  the  groove 


depth  we  assume  it  to  be  V-shaped  with  a  small  vertex  angle  6.  Then 

g (x)  =  N  dfc2  tan  g  (12) 

e.)  The  treatment  of  the  surface  diffusion  along  the  walls  of  grooves 
bears  a  close  resemblance  to  the  procedure  of  the  earlier  sections  dealing  with 
diffusion  in  the  grain  boundaries.  We  use  the  results  of  section  D  for  dg  to 
determine  Fsv,  the  quantity  analogous  to  Fv  of  section  (a) 

F  ---DA  2d  sec  g.  (13) 

sv  0  s  so  gov 


Here  Dg  is  the  surface  diffusivity  which  characterizes  the  free  surface.  (Our 
model  presumes  that,  without  grooving,  there  is  no  free  surface  because  of 
oxydizing  or  passivating  films.)  The  Aga  is  the  driving  force  coefficient  for  the 


surface  and  equals  Ag  cos  <j>0.  The  dg0v  is  the  groove  depth  for  the  o  segment  at 
the  V  vertex.  The  diffusion  contributions  to  (11)  have  been  omitted  since  the 
surfaces  of  the  groove  walls  are  smooth  and  there  is  no  preferred  direction  for 
matter  diffusion.  It  is  implied  that  there  is  no  electromigration-driven  pile-up 
of  vacancies  along  the  groove  wall  surfaces. 

As  before  the  Fsv  is  divided  into  fS(J,  which  are  now  all  equal  since  Ds 
is  independent  of  o.  The  back  distribution  into  the  segments  now  gives  rise  to  two 
different  possibilites :  segments  that  show  deficits  at  both  ends  and  those  which 
show  deficit  at  one  end  only.  Without  the  pure  diffusion  term  the  mass  motion 
equation  becomes 


3G  (x,t)  3G  (x,t) 

—?7 -  =  -D  A  — - 

3t  s  so  3x 


(2c) 


The  quantity  G(X,t)  is  that  groove  volume,  over  and  above  g(x,t)  previously 
defined,  which  is  acquired  by  virtue  of  the  surface  mass  motion. 

A  useful  long  time  solution  for  this  equation  is 

Gg ( x , t )  *  Bi  -  B2X  +  B2DsAst  (14) 

where  Bj  and  B2  coefficients  are  fitted  to  flux  conditions  at  the  vertices  or  zero 
flux  condition  for  those  grooves  that  terminate  in  a  segment  (single-ended 
grooves).  If  we  hold  to  the  V-shaped  groove  model  then  G  =  dg  tan  S'  where  now 
B'»S  and  S'  may  be  a  large  angle.  The  details  of  the  solution  of  this  surface 
aspect  have  not  yet  been  completely  worked  out. 


f.)  The  regions  of  material  excess  require  less  detail  than 
those  showing  deficit.  We  expect  to  select  a  criterion  for  the  breaking  of  the 
passivating  layer  by  internal  pressure.  After  this  is  exceeded  the  local  mass 
build-ups  with  time  will  cease  as  the  extra  material  is  extruded. 

F.  Present  situation  and  future  directions. 

Programming  has  proceeded  to  the  point  of  stage  c,  as  detailed  in 
section  E.  Trial  runs  on  a  single  grain  have  proved  satisfactory.  Further  work 

should  incorporate  the  remaining  stages  and  enlarge  the  operating  field  to  a  full 
stripe.  As  one  proceeds  with  ensuing  time  steps  the  grooves  become  incised  deeper 
and  surface  diffusion  replaces  grain  boundary  diffusion  as  the  principal  mode  of 
mass  transport.  For  the  cases  where  voids  have  developed  a  somewhat  different 
program  of  operations  will  be  needed.  It  is  quite  problematical  how  to  bring 
consideration  of  electrical  and  thermal  effects  into  play  but  they  certainly 


should  be  considered. 
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